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Abstract 

In this paper we present a new series of 3-dimensional integrable lattice mod- 
els with iV colors. The case N = 2 generalizes the elliptic model of paper [8]. 
The weight functions of the models satisfy modified tetrahedron equations 
with N states and give a commuting family of two-layer transfer-matrices. 
The dependence on the spectral parameters corresponds to the static limit of 
the modified tetrahedron equations and weights are parameterized in terms 
of elliptic functions. The models contain two free parameters: elliptic modu- 
lus and additional parameter rj. Also we briefly discuss symmetry properties 
of weight functions of the models. 



1 



1. Introduction. 



The purpose of this work is an enlarging of still poor zoo of 3D integrable 
lattice models. We say the model to be integrable if it possesses a family of 
commuting transfer matrices. The existence of such family is ensured, for 
example, by the construction of solutions of the tetrahedron equation, which 
is a three dimensional generalization of the Yang - Baxter equation [1-5] . As 
an example we can mention the N - color trigonometric model by Bazhanov - 
Baxter [6] . Although its integrability has been proved by a different method, 
the Boltzmann weights of this model are the solutions of the tetrahedron 
equation [7]. This solution, as well as the first known Zamolodchikov's one 
[1,2] (which is a particular case of the Bazhanov-Baxter model when N = 
2) can be parameterized in terms of trigonometric functions depending on 
tetrahedron angles. In our previous paper [8] we have constructed an elliptic 
two - color solution of a modified system of the tetrahedron equations, which 
provides the commutativity of two - layer transfer-matrices. This solution 
in a sense is not "full" and corresponds to the case of the so called "static" 
limit of tetrahedron equations. 

In difference from Bazhanov - Baxter model, for which the solution of 
tetrahedron equations contains six angular variables (five of them are in- 
dependent), the static solution of the modified tetrahedron equations from 
[8] can be parameterized by three angle - like variables and one additional 
parameter (the modulus of elliptic functions). 

In this paper we generalize this elliptic solution to the case of an arbitrary 
number N of spin values and obtain one more parameter, on which weight 
functions depend. This additional parameter is the same for all weights and 
new solutions are still static. Boltzmann weights of the new models like 
the weight functions of Bazhanov - Baxter model have the so called Body - 
Centered - Cube (BCC) form, invented by Baxter in [9] . This fact allows us 
to use the technique of the Star - Square relations, developed in [7]. 

The paper is organized as follows. In Section 2 we recall main definitions 
and notations, give the form of the Boltzmann weights and write out the 
modified tetrahedron equations. In Section 3 we formulate conditions for the 
BCC ansatz for weight functions to obey the modified tetrahedron equations. 
In Section 4 a natural parameterization of the obtained solution is given in 
terms of elliptic functions. In Section 5 we discuss symmetry properties 
of weight functions of the model. At last, Appendix contains a detailed 



2 



consideration of N = 2 case and an explanation of the transition to the 
weights of the model from paper [8]. 



2. Body-Centered-Cube (BCC) ansatz 
for weight functions 

In this section we recall some definitions from [7, 10] and give the explicit 
form of weight functions. We also write out modified tetrahedron equations, 
which provide a commutativity of two-layer transfer matrices [8]. 

Consider a simple cubic lattice C consisting of two types of elementary 
cubes alternating in checkerboard order in all directions (see Fig. 1). 
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Fig. 1 

At each site of C place a spin variable taking its values in Z N , for any 
integer N > 2 (elements of Zn are given by iV distinct numbers 0, 1, ... , N—l 
considered modulo N). To each "white" cube we assign a weight function 
W(a\efg\bcd\h) (see Fig. 2) 
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= W(a\eJ,g\b,c,d\h) 



Fig. 2 

and to each "dashed" cube - a weight function W(a\efg\bcd\h), where as 
usual a, . . . ,h - spin variables placed in the vertices of each elementary cube. 
Then the partition function of the model reads 

Z=J2 II W(a\efg\bcd\h) J] W(a\efg\bcd\h). (2.1) 

Spins '''white''' '''dashed''' 
cubes cubes 

Following [8] suppose that weight functions W and W satisfy the following 
equations: 

53 W(a A \c 2 , ci, c 3 |6i, 6 3 , b 2 \d)W'(c 1 \b 2 , a 3 , &i |c 4 , d, c 6 \b 4 ) 
d 

x W"(bi\d, c A ,c 3 \a 2 ,b 3 ,b 4 \c 5 )W {d\b 2 ,b A ,b 3 \c 5 ,c 2 ,c 6 \ai) 

= 53 W'"(bi\ci,c A ,c 3 \a 2 ,a4,a 3 \d)W (ci\b 2 , a 3 , a 4 \d, c 2 , c 6 |ai) 
d 

x W\ai\c 2) d, c 3 \a 2 , b 3 , ai\c5)W(d\ai, a 3 , a 2 \c±, c 5 , c 6 |6 4 ), (2.2) 

where W, W, W", W" &ndW,W ,W ,W are four independent pairs of 
weight functions. Suppose that the dual variant of (2.2) is also valid (with 
all VF's replaced by W's and vice versa). We call this pair of equations as 
a system of modified tetrahedron equations. Note that if put W = W in 
(2.2), we come to the standard version of tetrahedron equations (Eqs. (2.2) 
of [3]). In the Bazhanov - Baxter model all weights entering in tetrahedron 
equations are parameterized in terms of six angles 0i,...,9q satisfying one 
quadrilateral constraint. The explicit dependence on spectral parameters can 
be exhibited as (W = W) 

w - w(e 2 , e u e 3 ), w - w(ir - e 6 , e u 7r-e 4 ), 

W" -> W(9 5 , ir-e 3 ,n- A ), W"> -> W(6 5 , 6 2 , 9 6 ). (2.3) 
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In the static limit three spectral parameters 61,62,63 in weight functions 
W(6 1} 6 2 , 63), W(6 1} 6 2 , 63) are constrained by relation 6 1 + 6 2 + 63 = tt. Then 
formulas (2.3) are transformed as follows: 

W -> W{6 2 , 6±,tt 6\ 62), W > W{6 2 + 6 5 ,6 1 ,n-6 1 -6 2 - 6 5 ), 

W" -> W(0 5 , 01 + 62, TT - 0! - 62 ~ 5 ), W" -> W(0 5 , ^2, TT - 2 - 5 ) (2.4) 

and the same for W weights. Note that in this case the quadrilateral con- 
straint between tetrahedron angles is satisfied automatically. 

As it is shown in [8] the system of modified tetrahedron equations pro- 
vides a commuting family of two-layer transfer-matrices. Let us construct a 
horizontal transfer-matrix T(W, W) from alternating weights W and W as in 
Fig. 1. Here we imply that ingoing indices of the transfer-matrix correspond 
to the spins of the lower layer, outgoing indices - to the spins of the upper 
one and the summation over the spins of the middle layer is performed. As 
usual partition function (2.1) can be rewritten as 

Z = Tr{T(W,W)]} M , (2.5) 

where 2M is a number of horizontal layers of the lattice and we imply periodic 
boundary conditions. When equation (2.2) and its dual variant are satisfied, 
we have the following commutativity relation: 

T{W,W)T{W',W) = T(W',W')T(W,W). (2.6) 

Now let us specify the explicit form of weight functions. To make it, 
denote 

u = exp(2iri/N), lu 1/2 = exp(iri/N). (2.7) 

Further, taking x, y, z to be complex parameters constrained by the Fermat 
equation 

x N +y N = z N (2>8) 

and I to be an element of Z N , define 

w(x,y,z\l) = *[[—!L^. (2.9) 

In addition, define the function with one more argument 

w(x, y, z\k, I) — w(x, y, z\k — k,l G Z N , (2-10) 
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where 

$(/) = w '('+^)/2. (2.11) 

Let us mention also two formulas for w functions, which are useful for calcu- 
lations: 

w(x, y, z\l + k) = w(x, y, z\k)w(xu h , y, z\l), (2-12) 

w(x, y, z\k, I) = uj m /w(z, tu 1/2 y, ux\l, k). (2.13) 

Now introduce the set of homogeneous variables x iy i = 1, . . . ,8 and x 13 , 
%24, x 58 , x 67 satisfying 

T N _ N _ N N — N N N _ N _ N N _ N _ N (r, a \ 

x 13 — x i 3 ' x 24 — x 2 x 4 i x 58 — x 5 x 8 > x 67 ~~ x 6 x 7 ■ v^' 1 V 

Using all these notations we define the weight function W (a\ef g\bcd\h) as 

IU a e,/,</ 6, c, d \h = ^ ~ ( , : — w — i , , , , v 

atz N w{x 8 ,x 58 ,x 5 \e,c + a)w{x 7 /u,XQ 7 ,x 6 \f,b + a) 

(2.15) 

In fact, the Boltzmann weights of form (2.15) generalize the weight func- 
tions of the Bazhanov -Baxter model [6, 10]. Up to inessential gauge factors 
the latter corresponds to the choice 

x 5 = xi, x 6 = x 2 , x 7 = x 3 , x 8 = x A . (2.16) 

Following paper [6] we will call formula (2.15) as a Body-Centered-Cube 
(BCC) ansatz for weight functions. 

3. Star-square relation and the proof of mod- 
ified tetrahedron equations 

In paper [7] tetrahedron equations for the Bazhanov - Baxter model with N 
states were proved using the so called Star-Square and "inversion" relations. 
We will follow the method of this paper for weight function (2.15). 
First recall "inversion" relation for functions w(x, y, z\k, I): 

v w(x,y,z\k,l) (l-z/x) 
k ^ N w(x,y,u;z\k,m) 0l > m (l- z»/x»)' ^ } 
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where I, m e ^at, x,y,z satisfy (2.8) and 5^ m is the Kronecker symbol on 
Zn- 

To write down the Star-Square relation introduce a non-cyclic analog of 
w function, defined recurrently as follows: 

w(x\0) = 1, . = jr-, Z G Z, (3.2) 

V 1 ; ' w(x|Z-l) (1-W)' ' V ; 

where Z is the set of all integers. It is obvious that 

w(x,y,z\l) = (-) l w(x/z\l), leZ N , (3.3) 

where index /, being considered modulo N, is interpreted as an element of 
Zn. Then we have the following identity 



w(xi, zi|q + a)w{x 2 , 2/2, ggjo + 
tz N w(x 3 , J/3, z 3 \c + a)w(x 4 , t/4, z 4 |rf + a) J 

(x 2 zi M ^2 ) 0-6 (xi 2/2 /ar 2 *i ) 6 (z 3 /ya ) c (z 4 M ) d 



$(a - 6)c>+ 6 )/ 2 

w(u;a;3a;4ZiZ2/xia;2-23-24|c + d - a - b) (li 

x /X'42i,, : ,x' 3 2 2 | n ^3^1 1 : / x '4^2,, 77 ' 
u>( 1 a — a)w{ |c — o)u;( \c — a)w{ \a — b) 

X\Z 4 2-1^3 X2Z4 

where the lower index "0" after the curly brackets indicates that the l.h.s. of 
(3.4) is normalized to unity at zero exterior spins, and the following constraint 
is imposed 

V\ 2/2 z 3 z 4 . 

= uj. (3.5) 

Z\ z 2 2/3 Va 

We call relation (3.4) as the Star-Square one. Note that the separate u>'s 
in the r.h.s. of (3.4) are not single-valued functions on Zn, while the whole 
expression is cyclic on the exterior spins a, b, c, d. 

The proof of (3.1), (3.4) was given in [7] and we refer the interested reader 
to this paper. 

Now we turn to relation (2.2). Instead of weights W, W, W" , W" let 
us substitute into (2.2) explicit formula (2.15) with corresponding sets of 
parameters: 

W, W, W" , W" - W(x h Xij ), W(x[, a/^ww, 4), W(x'l', 4), 
W,W',W",W" - Wfaxij), W&rtj), Wix';,^), W (3.6) 
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Let us multiply both sides of (2.2) by the following product of w weights: 





c 4 ,a 2 + l 2 ) w(x' 7 ,ujx' 67 ,x' 6 


a3,C 4 + l 3 ) 


1 — 111 1 — /// — 

w{x 3 /u, x 13 , X 


111 

1 


C 6 , ay + ly) W(X4, X 24 , UX 2 


a 4 , c 3 + h) 



X 



X 



/ — II — // — // / 

w{x 8 ,x 58 ,x 5 /uj 


b 2 ,c 2 + m 2 ) w(ujx 8 , x' 58 , x' 5 \ 


c 2 , b 3 + m 3 ) 




by,C 3 + my) W{X 3 /CV, Xy 3 , Xy\c 6 ,b 4 + m 4 ) 



(3.7) 



and sum over ay, a 2 , a 3 , 04, by, b 2 , b 3 , 64. Note that due to "inversion" relation 
(3.1) we do not lose any information after such transition from spins a^, bi to 

The functions w in expression (3.7) are chosen in such a way that using 
relation (3.1) we can calculate the sums over the spins ay, a 2 , a 3 , 04 in the l.h.s. 
and those over by,b 2 ,b 3 , 64 in the r.h.s. of the obtained equation and cancel 
the summations over the spins cr's, which come from expression (2.15) for the 
functions W's and Ws. Now let us consider applicability conditions (3.5) of 
Star-Square relation (3.4) for the sums over ay, a 2 , a 3 , 04 in the r.h.s. and over 
by,b 2 ,b 3 , 64 in the l.h.s. of the obtained equation. We obtain eight conditions 
on re's and rr's. Applying relation (3.4) eight times and calculating the sums 
over d spin in the l.h.s. and r.h.s (the spin structure of the sums over d has 
the form of relation (3.1) and we demand that corresponding variables x, y, z 
entering in the arguments of functions w(x, y, z\k, I) are constrained in such 
a way that relation (3.1) can be used), we come to the equation without any 
summation. The l.h.s. and r.h.s. of this equation consist of the products of 
ten w functions with the same spin structure and expressions like x™ 2 coming 
from relation (3.4). Let us impose all necessary constraints on parameters 
Xi, Xij, Ti and Tij to satisfy this equation. On obeying of these constraints 
one can show that spin independent multipliers coming from relations (3.1), 
(3.4) coincide. 

We also have to satisfy the dual variant of (2.2). Hence, we must add 
a dual set with x replaced by x and vice versa to all obtained constraints 
on parameters x and x . We will not write out here all these relations. A 
detailed analysis shows that we have two solutions. The first one corresponds 
to the choice 

Xi = Xi, Xy 3 = Xy 3 , X 24 = X 24 , X 58 = X 58 , X 67 = X 67 , 1 — 1, . . . , 8 (3.8) 

and the same for x', x" , x'" . This choice corresponds to the Bazhanov - 
Baxter model, considered in [7]. 
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But there is also another possibility. It will be convenient to fix a nor- 
malization of all parameters x in W functions as 

x 3 = 1, X4 = 1, X7 — 1, x$ = 1 (3.9) 

for all sets x, x' ', x", x'" and x, x', x", x w . Then all parameters x can be 
expressed in terms of x as: 

Xi = l/x 2 , x 2 = X 5 = l/x 6 , X 6 = l/x 5 , 

X 13 = C^fM = w l/2^ > ^ 5g = w l/2^ > ^ = W -l/2^. (3.10) 

Let us introduce the following notations: 

_ X1X2 f _ X 5 , X_6 y. X 58 X 6 7 

X5X6' 1 WXi' t2 Xi ' ^ X13X24' 

_ x % ^13^58 a _ , ,^5 Xi 3 X 67 _ ~ ~ (q 1 1 \ 

~ XTX 67 X 2 4' 32 ~ U WX 58 X24> J3-X 5 X 6 . [6.11) 

and the same for the sets of Then all constraints on parameters 

take the form 

i _ 4-in + — +1 + — +11 +1 — +111 +1 — +11 +11 — +111 

L l — L 2 1 L 2 — t 2 , ^3 — t 2 > l \ — ^3 1 L 3 — L 3i l \ — l \ 1 

3i= h=f 2 , J3=J^, j[=fi', j' 3 =%, Ji=Ji', (3.12) 

Q rJM /->-» „ iy*^ iyf nfJ^I 

■j — o — o — o j ^24 "^24 — 24 24 " 

Relations (3.9-3.12) together with consistency relations (2.14) are sufficient 
conditions for weight functions (2.15) to satisfy modified tetrahedron equa- 
tions (2.2) and their dual variant. In the next section we will obtain a natural 
parameterization for relations (2.14), (3.11-3.12) in terms of elliptic functions. 



4. Parameterization 

Recall that consistency relations (2.14) connect iV-th powers of Xj and x^-. 
Hence, it is convenient to consider iV-th powers of (3.11-3.12) and introduce 
new parameters 

S 2 = s N , T? = -t»/S, J? =jf, i = l,...,3. (4.1) 

From (3.11-3.12) we see that parameter S is the same for all weights entering 
in relation (2.2). 
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Parameters S, Tj, J«, i = 1, 2, 3 are defined by four independent variables 
xi, %2, £5, xq. Therefore, there are three additional constraints between these 
values. After simple calculations we obtain 

, = Tih-T 2 J 2 {Tt-Tpjl + SJMT,) 

J?> TxJa-TaJx' 3 (1 — T 1 2 T 2 2 )(7i^2 — T 2 Ji) 1 ] 

and 

fs g± jf I M ± gfl i i = 12i (43) 

where we introduce a new parameter P. 

Note that definitions (4.1) contain only the second powers of variables S, 
T{, Ji and in (4.2-4.3) we have chosen some signs in a convenient way. Using 
(4.2) it is easy to obtain that we can add the case i — 3 to relations (4.3). 
The validity of formulas (4.2-4.3) can be checked by substitution of relations 
(3.11), (4.1) and consistency conditions (2.14). 

Introduce an elliptic curve 

p = S + x 2 -Sy 2 (l + Sx 2 ) ^ 
xy 

Points (Tj, Jj), i = 1, 2, 3 belong to this curve. We can uniform it in terms of 
elliptic functions (see, for example, [11]). Note that formulas (4.2) take the 
form of addition theorems for elliptic functions. Then we obtain 

T^imWmM, Jj = Sa{ ^ i \ j = 1, • • • , 3, 

P = — 2im 1 / 2 sn(r))cn(r))dn(i]), S = msn 2 (i]). (4.5) 
where sn, cn, dn are elliptic functions of modulus m and 

u 3 = u l +u 2 . (4.6) 

In such a way we can obtain the parameterization for xf and xfy When we 
take roots of iV-th powers, a phase ambiguity appears. Making an appropri- 
ate choice of these phases so that (3.11-3.12) to valid we obtain the following 
formulas 



Xl = uo- 1 ' 2 



sn(?7 — Uz)sYi{rf) ~]}^ N 
sn(wi)sn(w 2 ) J 
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x 2 = c// 2 {m 2 sn(?7 — ■u 3 )sn(ry)sn(Mi)sn(M 2 )} 1//Ar , 

x _ ^,i/2( sn(txi)sn(7y - m 3 ) ) 1/jV ^ = „-i/ 2 ( sn(^ 2 )sn(?7 - u 3 ) 1 A 
5 \ sn(-u 2 )sn(?]) J ' '' " \ sn(ui)sn(?7) 



(4.7) 



and 



x 13 = UJ 



-1/2 f e(u 3 )e(0)H(rj - Ul )H(rj - u 2 ) 
1 8(77 - u 3 )e(rj)H( Ul )H(u 2 ) 



l/N 



_ 1/2 r 9(m 3 ) 9(0)9 (77 - ^1)6(77 - ug) ~i i/iv 
X24_a; 1 9(77- M3 )9(77)9( Ml )9(w 2 ) ' ' 

1 9(77 -u 3 )H(rj)e( Ul )H(u 2 ) 



rr 58 = cj 1/2 



_ 1/2 f jJ(^ 3 )9(0)9(77 - (77 - tx 2 



l/JV 

X67_ I 9(7/- M3 )iJ(77)i/( Ml )9( M2 ) J ' ^ 

where H{u) and Q(u) are Jacobi elliptic functions (see [11]). Suppose that 
real parameters ui, u 2 , 77 satisfy 

< ui.2,3 < 77 < (4.9) 

where Ai is the complete integral of the first kind of the modulus m. These 
conditions guarantee that all values of the elliptic functions in (4.7-4.8) are 
non-negative, and we choose the positive values of all roots of iV-th power in 
formulas (4.7-4.8). 

Now we can rewrite the Boltzmann weight as a function of four new 
parameters: W(ui, u 2 \m, 77) (we omit the dependence from spin variables). 
It is on easy to check that dual weights W can be obtained by the shift of 
parameter 77: 

W(ui, u 2 \m, 77) = W(ui, u 2 \m, 77 + iM'), (4-10) 

where Ai' is the complete integral of the first kind of the complementary 
modulus m' = y/l — m 2 . Equations (3.12) relate the sets of new parameters 
in different weights. Then we obtain 

W^W(u 2 , Ul ), W ^W(u 2 + u 5 , Ul ), 
W" -> W(u 5 , u x + u 2 ), W" ->W(u 5 ,u 2 ). (4.11) 

Parameters m and 77 are the same for all weight functions and we omit them 
in (4.11). Also note that we have chosen three independent parameters in 



11 



(4.11) as ui, u 2 and u§ to emphasize the connection of our parameterization 
with static limit one (2.4). 



5. Symmetry Properties 

In this section we discuss symmetry properties of weight functions (2.15) 
of the model with respect to the group G of transformations of a three- 
dimensional cube (see [10, 12]). Rcall briefly some definitions. Group G has 
two generating elements p and r. Any element a of G can be expressed as a 
composition of these two elements. We define the action of elements r and 
p on the set of spins {a\e, f, g\b, c, d\h} as follows 

r{a\e, f, g\b, c, d\h} = {a\f, e, g\c, b, d\h} (5.1) 

and 

p{a\e, f, g\b, c, d\h} = {g\c, a, b\f, h, e, \d}. (5.2) 

Further it will be convenient to remove normalization conditions (3.9) and 
to restore the homogeneity of the parameters re's. 

It is easy to see that weight function (2.15) is invariant under the action 
of element r, if we make the following transformation of parameters: 

(x 3 , X 13 , Xx) -> (x 3 , X 13 , Xi), (x 4 , x 24 , x 2 ) -> (x 4 , x 24 , x 2 ), 
(x s , x 58 , x 5 ) -> (x 7 /u, x 67 , x 6 ), (x 7 , x 67 , x 6 ) -> (ux s , x 58 , x 5 ). (5.3) 

The action of element p is less trivial and can be obtained with the help of the 
Fourier transformation (see [10]). Introduce two more additional parameters 
u and v which are defined as 

N N N N N N N N 

rv* ty ty ty ry * ry^ ' ry * ry A * 

N _ x 3 x 5 x l x 7 N _ x 2 x 8 x 4 x 6 /r A \ 

a — N , v — N . {o. 4) 

A phases of parameters u and v can be chosen in an arbitrary way. Then 
one can show that the following identity is valid: 



r ^ w(x 3 ,xi 3 ,xi\d,h + a)w(x4,x 2 4,x 2 \a,g + a) 
[ a ^ N w(x 8 , x 58 , x 5 \e, c + a)w(x 7 /u, x 67 , x 6 \f,b + a) 



_ , .dh+ch-bf-bg w(x 4 X6, x 2 v, x 2 , x 8 \a + b, f + g) ( , 



w(xix 7 , X\U, x 3 x 5 \e + h,d + c) 
r ^ w(x 7 x 13: xiu, x 3 x 57 \e, d + a)w(x 2 x 68 , x 2 v, x 6 x 24 \g, b + a) 
l aez N w^x^.ujxiu.uxxx^c, h + a)w(x 4 X68, x 2 v, x 8 x 2i \a, f + a) 
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where the lower index "0" after the curly brackets implies that the expres- 
sion in the curly brackets is divided by itself with all exterior spin variables 
equated to zero. 

The arrangement of spins in the sum of the r.h.s. in (5.5) corresponds 
to the p-transformed set of the indices (see (5.2)). Hence, we can define the 
action of p transformation on the parameters as 

(X 3 , X 13 , Xi) A (x 7 Xi 3 , X ± U, X 3 X 57 ), (X 4 , X 24 , X 2 ) A (x 2 X e8 , X 2 V, X e X 24 ), 
(x 7 ,X 57 ,X 5 ) A (x 5 X 13 ,UJXiU,UJXiX 57 ), 

(x 8 , x 68 , x%) A (UX4XQS, x 2 v, x 8 x 24 ). (5.6) 

Using (5.5) we can choose gauge factors before formula (2.15) in such a way 
that the whole W function will be invariant under the transformations from 
the group G. We will not write this symmetric form of weight functions here. 
Note only that after p transformation (5.6) W function is transformed into 
W and vice versa. Using relations (4.7-4.8) from the previous section we can 
write the following transformation rules for parameters Ui,u 2 , r) with respect 
to the elements r and p from the group G: 

Ul A U 2 , U 2 Auj, 7] A 77, (5.7) 

Ui-^iM'-Ui, u 2 A m + u 2 - iM', r]-^r] + iM f . (5.8) 

We hope to use these symmetries for the calculation of the partition function 
for our models. 

Appendix. A detailed consideration for the 

case N = 2. 

The BCC form of the weight is useful for proving the modified tetrahedron 
equations. But for other purposes a symmetrical form of the Boltzmann 
weights, which differs from the BCC one by some gauge transformation, is 
more preferable. In this section we find this gauge multipliers for the case 
when N = 2 and obtain a table of Boltzmann weights analogous to Baxter's 
table from [9] and endowed by obvious symmetrical properties. Also we show 
that fixing r\ = M. we obtain the particular case considered in [8]. 
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We start from a simple gauge transformation of our BCC weight (2.15) 
and consider the weight W: 



W'(a\e, f,g\b,C,d\h) = ^-ifi+ec+a+h^af+ch-dh+ag+cg-df x 

x Yf' d ' h ' c Y?^ c : h ' h ' 9 W{a\e, f, g\b, c, d\h), 



Y nfhn Y f 



a,.f>b,g 1 b,h,djY e ^j^ a 



where 



.7T 



Y a ,f,b,g = exp{-i-fg}exp{mab(f - g)} 



(A.l) 



(A.2) 



It is convenient to introduce spins (— l) a = ±1 instead of a = 0, 1. Further 
we shall use only the multiplicative spins and write a instead of (— l) a , etc. 

Absolute values of these weights ( both W and W) depend only on three 
Baxter's sign variables abeh, acfh, adgh, abed. Therefore up to signs there 
exist sixteen different weights (A.l). In terms of multiplicative spins we can 
write the following parameterization: 



W' = £W s x 
x exp{(Mi + N x )agbf + (M x - N 1 )cedh+ 
+ (M 2 + N 2 )agce + (M 2 - N 2 )bfdh+ 
+ (M 3 + N 3 )aefg + (M 3 - N 3 )cgbf+ 
+ (M + N )abcd + (M - N )efgh}, 



(A.3) 



where, calculating all sign factors, we parameterize Ws as follows: 



abeh 


acfh 


adgh 


W s (a\efg\bcd\h) 


+ 


+ 


+ 


Po 




+ 


+ 


Ri 


+ 




+ 


Ri 


+ 


+ 




Rz 


+ 






abPi 




+ 




acP 2 






+ 


adP 3 








abcdSo 
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Parameter £ in (A. 3) is a normalization factor and we choose it so that 
Po = 1- Now we can obtain exp(16Mj), exp(16Aj) and W§ just by multiply- 
ing and dividing BCC weights for several sets of the spins. 

It is useful to express the answer in terms of elliptic functions of another 
modulus k instead of the used functions of the modulus m 

k = . (A.4) 

1 + m v ' 

Denoting the complete elliptic integrals of the first kind of the moduli m and 
mf = yl — m 2 as M. and Ai', and of the moduli k and k' = yl — k 2 as K. 
and K,', we have the connection between them 

K! = (l + m)M, 2/C = (l + m)M'. (A.5) 
One can easily obtain the reparameterization of x 2 and x?- using 

?m 1 ^ 2 sn(t>, m) = fc' — S / ' 7 — — , where — - — = iv. (A. 6) 

cn(x, k)dn(x, k) 1+m 

To calculate expressions for W S) exp(2Mj) and exp(2Aj), we need three more 
complicated identities: 

sn(w 1 + v 2 ) - sn(vi) - sn(v 2 ) - msn(v 1 )sn(v 2 )sn(v 1 + v 2 ) 
sn(vi + v 2 ) + sn(vi) + sn(v 2 ) - msn(-ui)sn(w 2 )sn(t;i + v 2 ) 

sn(xi)sn(x 2 ) 

cd(x!)cd(x 2 ) ' 

sn(wi) + sn(t> 2 ) + sn(f 1 + v 2 ) + msYi{vi)sYi{v 2 )sYi{vi + v 2 ) 
sn(f 1) + sn(t> 2 ) + sn(t>! + v 2 ) — msn(f 1 )sn(v 2 )sn(w 1 + v 2 ) 

cd(xi + x 2 ) 

cd(xi)cd(x 2 ) ' 

sn(i>i) + sn(f 2 ) — sn(t>i + v 2 ) — msn(vi)sn(v 2 )sn(t;i + v 2 ) 
s\i{vi) + sn(f 2 ) + sn(f x + v 2 ) — msn(v 1 )sn(v 2 )sn(u 1 + v 2 ) 
= — /c 2 sn(xi)sn(x 2 )cd(xi +x 2 ), 

(A.7) 

where in the LHS we imply the modulus m, and in the RHS - k, the argu- 
ments Vi and Xi being connected by (A. 6). 
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We would like to obtain a model with real weights Ws- To do this we 
have to consider a regime when ui, u 2 , w 3 are pure imaginary satisfying 
< Im« li2)3 < M' and rj = M. + it, a small e being real. Defining further 
the variables -21,2,3, A by 



2^1,2 



= -wi, 2 , z 1 + z 2 + z 3 = K, 



2A 



1 + m ' " 1 + m 

and using repeatedly (A. 7), we obtain the exponentials 

cxp 2JVi = e ^ /4 {-iA;sn(A + zi)cd(A + 2i)} 1/4 , 

cxp 2AT 2 = {-iksn(X + z 2 )cd(\ + z 2 )} 1/4 , 

exp2N 3 = {-iksn(\ - z 3 )cd(A - ^ 3 )} 1/4 , 

exp2A" = {-iA;sn(A)cd(A)} 1 / 4 ; 

r - c d(A + zi) 1 1/4 , .cd(A + 2 2 ) , 1/4 

exp 2Mi = {1 77— 7} > ex P 2M 2 = {2 } 1/4 , 
sn(A + Zi) sn(A + z 2 ) 



(A.8) 



exp2M 3 = {. Cd ^ Z3 \ }^, exp2M = {- 
sn(A - z 3 ) 



.sn(A) 1/4 
'cd(A) 1 ' 



(A.9) 



(A.10) 



In this formulae Zi are real, < Zi < K, and A = i/C'/2 + t' , t' being small 
and real. The expressions contained in the curly brackets belong to the first 
or fourth quadrants of the complex plane, the fourth power roots of these 
values are defined so that they belong to the same quadrants. As to Ws, it 
coincides with the elliptic solution of ref. [8]: 



Po = l, Pi 



sn(zj)sn(z k ) 
\| cd(zj)cd(z k y 



I SI1 ( Z ' J 

S = Aysn(, lM ,2)sn(, 3 ), Ri = ^ cd{z ^ {zk) , (A.ll) 

where i,j,k are any permutation of 1,2,3. Here the square roots of the 
positive values are supposed to be also positive. Our weights coincide with 
those of [8] when 

A = iy ^r] = M. (A.12) 

The dual weights W can be obtained from these ones by the changing 
Mi — > —Mi, i = 0, ...,3, and So — > —So, so that the terms Mi, 2)3 and A^i j2 , 3 
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are the pure gauge of the system of the modified tetrahedron equations and 
the only difference with the case (A. 12) consists in the Mq and Nq terms. 

Note that the extraction of the A - dependent terms exp(M ) and exp(iV ) 
as common multipliers of the A - independent weight Wg is the property of 
the choice N = 2, when N ^2 this property fails. 
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